Fetching marked items from an unsorted database in NMR ensemble computing 
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Searching a marked item or several marked items from an unsorted database is a very difncult 
mathematical problem. Using classical computer, it requires 0(N = 2 n ) steps to find the target. 
Using a quantum computer, Grover's algorithm uses 0(V N = 2 n ) steps. In NMR ensemble comput- 
ing, Brüshweiler's algorithm uses log N steps. In this Letter, we propose an algorithm that fetches 
marked items in an unsorted database directly. It requires only a single query. It can find a single 
marked item or múltiple number of items. 



Unsorted database search is a very important mathematical problem and it has many applications. The problem 
can be expressed as follows: Suppose we have a query function f. f is known to be zero for all inputs x, /(x) = 0, 
except for one or some of the items in a database containing N — 2™ items, say, x = z and f(z) — 1. The problem is 
to find z. On a classical computer one has to evaluate the function f about 0{2 n ) times. Using Grover's algorithm 
0, we can find z in 0(2™/ 2 ) steps which is optimal for pure or effective-pure states [||. In contrast, Brüschweiler's 
algorithm finds z in 0(n) steps [|| in an NMR ensemble computer by using truly mixed spin states. The speedup 
I , is achieved by the massive parallelism of representing input states using molecular sub-ensembles with diffcrcnt spin 
' states. Different spin states perform diffcrcnt computations simultaneously in molècules. This parallelism is classical 
ÇS| , in nature 0. In this Letter, we present an algorithm that uses the topological nature of the NMR spectrum to 
• actually fetch the marked items from an unsorted database. Like the Brüschweiler algorithm, our algorithm is also 
for ensemble computers. The algorithm is a hybrid algorithm that combines both DNA computing and quantum 
computing. 

Spin density operator is used to describe NMR ensemble system Js],^). The usage of polarization operators permits 
a concise and suggestive description of spin ensemble quantum computation. Spin systems containing n spin 1=1/2 
nuclear spins are chosen to encode the database. The eigen-state of the nuclear spin systems \4> m ) of a weakly-coupled 
system can be expressed as a direct product of the single-spin eigen-functions o^}, \f3i) of Ii Z , for instance, 
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\(f> m ) = \a 1 )\a 2 )\0s).~\a<n-i}\Pn) = \aaf3...af3) = |001...01>, (m = 1, ...,2 n ), (1) 
wherc \a) = |0) and |/3) = |1). In Liouville space, the 2™ pure states can be expressed by direct products of the 

I a = \a){a\ = \{l + 2I Z )=(\ ° Y (2) 



i = i/w = 5(1-2/,) 



(3) 



The corresponding relationship between a density operator p and the pure state \<j>) is 

\4> m ) = \aa(3...a(3) p m = / 2 "/f ...7«_ x /f (4) 

Before presenting our algorithm, let's first briefly review the general description of a Fourier spectroscopy Sup- 
pose before the acquisition, the density matrix of a spin system is p(0-). In order to acquire the wanted information, 
a proper radio frequency (r.f.) pulse such as non-selective r. f. pulse is operated on p(0_) 

p(0+) = Ry(fif)p(0-)Rv(P) (5) 
where R y is the rotation operator along y axis. 

R y {[3) = exp{-ií3F y ] (6) 

where the F v = Y] Ik y ■ The subsequent free evolution is 

Í^. = -i[H,p]-t{p(t)-p } (7) 
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where H is Hamiltonian and T is relax super-operator. With the Liouville super-operator L = —iH — T and the 
equilibrium dcnsity operator po commuting with the unpcrturbcd Hamiltonian H, The above eq. (7) can be rewritten 
as 

^ = L{ P (t) Po } (8) 

The solution of this equation is 

p(t) = exp{Lt){p{Q+) - po} + Po (9) 
The magnetization M + (t) is proportional to the value of the operator F + = F x + iF y = Ikx + i Iky 

M+(t) = Njh(F+)(t) = Njhtr{F + p(t)} (10) 

where N is the numbcr spin systems per unit volume. 

In a high magnètic field, the Liouville operator L is invariant under rotation about the z-axis and does not mix 
components of p belonging to diffcrcnt coherence orders p, then a general representation of the complex free induction 
decay is 

M+(t) = N 1 htr{F+exp{Lt)p{Q + )} (11) 

In the absence of degenerate transitions, each off-diagonal matrix element of p(t) in the eigenbase of H evolves 
independently 

M+{t) = N 1 nY J F+Prs^+)exp{{-iw rs - X rs )t} = N-fhJ2 F+p sr {0 + )exp{(iu rs - X rs )t} (12) 

rs rs 

where the transition freqüències are 

u) rs = H rr - H ss = (r\H\r) - (s\H\s) (13) 

and the relaxation rates are 

-Vs = —Fl rs rs = (14) 

The complex Fourier transformation of the frec induction signal yields the complex spectrum 

S(w) - JV 7 ft£ F "M0 + ) . A 1 (15) 

rs 

Obviously, from this formula of complex spectrum, peaks of the spectra have an one to one correspondence to the 
transitions between energy levels. The states of an ensemble computer can be identificà by transitions between energy 
levels in contrast to the common correspondence between the ensemble- computer states and the energy levels themselves. 
As we will show later that peaks of the spectra can be used to label the states of an ensemble computer. A main theme 
of our algorithm is the use of this correspondence. 

Firstly, wc will elucidate the correspondence. The Hamiltonian in a weakly-coupled system can be expressed as 

H = ^{ujl jz + ^Jjkljzhz) (j = 1, 2, ....n) (16) 

3 k)j 

The eigen-state |r) can be expressed as a direct product of the single-spin eigen-functions 

\r) = \hÍ2-i n ) ii^Oorí (17) 

The eigen- value is 



(r\H\r) = (iiÍ2..A n \Y(u j I jz +^2nJ jk I j J kz )\i 1 i 2 ...i n ) = ( (-lfiuj + '£*J j k(-l) ii+ik 

3 k)j 3 V k)j j 



(18) 



we adapt the convention \ii) = |l/2) = |0) corresponding to excited state (high frequeney) and \ii) = \ — 1/2) = |1) 
corresponding to ground state (low frequeney). Generally speaking, in one dimension spectrum, only singlc-quantum 
coherence can be observed directly. The transition selection rules are AM = where M is magnètic quantum 
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number. In our algorithm, we use an aucilla bit which is labled as zeroth qubit: Iq. The energies of the system when 
the aucilla bit is at and 1 are, 

(0,i 1 ...i n \H\0,i 1 ...t n ) = ^{ujo + Y / (- 1 Y J ^+Y, 7rJ ^- i y k +E E KJjk(-l) ií+Íh } (19) 

j=l fe=l 3 k)j=l 



(l,z 1 ...i„|ií|l, il ...z„}-Í[- Wo + ^(-l) í ^,-^7rJ ofe (-l)^+^ Yl «M-lY í+ % (20) 

j=l fc=l 3 k)j=l 

respectively. If we observe the spectrum of the aucilla bit Iq , the transition freqüències will be 

cj + J2^Jok(-l) ik , (fc=l,...,n). (21) 
fc=i 

In the spectrum there are 2™ peaks. Each peak corresponds to a transition between \0ixÍ2···in) and \lixÍ2···in) in the 
aucilla bit's spectrum. The frequeney of the peak is determined by eq.(pï|). If the state of the aucilla bit is in 0, then 
the transition is from \0ixÍ2···i n ) to \liiÍ2-..i n ), and the peak is upward in the spectrum. If the state of aucilla bit at 
the aequisition is at 1, then the transition is from \lixÍ2···in) to \0ixÍ2···in), and the peak is downward. Thus the 2™ 
number of peaks in the aucilla bit spectrum correspond to 2" numbers: \iiÍ2-..i n ). The upward or downward nature 
of the peak indicate the aucilla bit's state before the aequisition is or 1. The \i\Í2---i n ) = |00...0) state is the far 
left peak(highest frequeney) in the spectrum, and \i\Í2---in) — 1 1 1 ... 1 ) is the far right peak in the spectrum (lowest 
frequeney). In between them are the other states of the system. 

With the above knowledge in mind, it is easy to present our algorithm: 1) first prepare the n bit with the aucilla 
bit system in the Iq state. This state actually represents 2" items of the database, since 

I§ = I%hh...In = W + + + %) 

= Wï-tf + + - + I l I 2- I n): (22) 

2) apply the query function / to the system and store the results of the query in the aucilla bit. The query function 
corresponds to an unitary transformation 

f = Tr{UIiï Pm U t I 0z }. (23) 

For truly mixed state, 

M M 

f = '£F(^P j )=F(J2 I oPJ^ (24) 

3=1 i=l 

3) measure the aucilla bit's spectrum and read out the results from the spectrum. Marked items will be those states 
with peaks downwards in the spectrum. 

We illustrate our idea using one simple example with n = 2. The unsorted database can be expressed as (the extra 
bit is ineluded in) {000 (/^f/f), 001 (iff/fif ), 010 (I Q /f/ 2 Q ), 011 (ifflfif )}• The energy levels and transitions 
of the extra bit are: (000|ií|000) = 2 (^o + uji + ui 2 + ttJqx + ttJq 2 + kJi 2 ), 

(100|Jf|100) = -|(o; + uji + uj 2 - ttJoi - 7rJ 02 + irJ 12 ), (001|iï|001) = |(w + wi - w 2 + ttJ i - 7rJ 02 - 7rJi 2 ), 
(101|íf|101) = -5(0)0 +U1-W2- ttJoi + ttJ 02 - 71-J12), (010|íf|010) = 2(^0 -W1+U32- ttJqi + 7rJ 02 ~ ttJi 2 ), 
(110|ií|110) = -3(<xJ - w x + uj 2 + ttJoi - ttJ 02 - nJ 12 ), (011|ií|011) = \{ujq - ujx - uj 2 - ttJ i - iïJq 2 + %J 12 ), 

(lll|i/|lll) = -5(0)0 - Ux - LÜ2 + TT^Ol + 7TJ 02 + 7rJl 2 ), CJ 4 = ^0 + 7t(J 1 + J 02 ), W].5 = ÜJq + 7t(J 1 - J 02 ), ^26 = 
Uq + 7T(- Joi + J02), ^37 = ^0 + 7t(-J 1 ~ Jq 2 )- 

There are four transitions which corresponding to {00(/f 01 (7f if), 10(/f/ 2 a ), jf )}. When we aequire 

spectrum on the state Iq, all the peaks in the spectrum are up. After performing the oracle query /, the peaks 
corresponding to marked states will be downwards. 

It will be easy to implement the above example in an NMR system. The physical system can be the C labeled 
alanine "C* 1 ^ - 13 C H(NH£) - 13 C 2 OOH. The solvent is D 2 0. The J couplings are: J i = 35.1Hz and 
Jq2 = 54.2ÜZ. Suppose the marked items are 10 and 11. Starting from thermal equilibrium state: 

a(0_)=4° + 4 1 +4 2 , (25) 

state Iq is prepared. It can be achieved by a sequence of selective pulses, non-selective pulse and J-coupling evolution: 

^Qrad (26) 
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In which, subscripts denote the directions of the radio frequency. Superscripts denote the nuclei which the radio 
frcqucncy is operated on. Grad refers to gradient field. Then the query function operation is performed. It is 
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(27) 



The corresponding pulse sequence is 



* y 



(-)" 

v o 'y 



(28) 



where r = tt4— . 

Wc acquirc the spcctrum of the aucilla bit after the U transformation is operated on the state Iq . The peaks down 
correspond to the marked states 10, 11 respectively. The left one down is corresponding to 10, and the right one 
down is 11. The spectrum from computer simulation is shown in Fig.l. 

In summary, we have proposed an algorithm that finds marked states in an unsorted database with a single query. 
It requires much less computing time and is thus more robust against decoherence. In addition, it has several other 
advantagcs. The algorithm does not require the prior knowlcdge of the number of marked states in the database. It 
uses the shape of the spectrum, which is a topological property of the spectrum to read out the marked states. This 
makes it more robust to errors and has much less demand on the accuracy. As for practical implementation, it is very 
important to choose the right working media. In particular, the aucilla bit should have J coupling with all the other 
bit. 
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Figure captions: 

Fig.l. The aucilla bit spectrum after performing query. The downward peaks correspond to the marked item 10 and 11. 
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